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Abstract 

Limits on the regions of and x where the operator product 
expansion can be safely used, at small values of x are given. For 
a fixed large there is an xo((5^) such that for Bjorken x- values 
below xq the operator product expansion breaks down with significant 
nonperturbative corrections occurring in the leading twist coefficient 
and anomlous dimension functions due to diffusion of gluons to small 
values of transverse momentum. 



1 Introduction 



For hard scattering processes involving two transverse momentum 
scales, one of which is the hard scattering scale, factorization and the 
Dokshitzer, Gribov, Lipatov, Altarelli, Parisi (DGLAP) equation [1-3] 
furnish the basis for a systematic description of the dependence of 
cross sections on the hard scale. For example, in deep inelastic lepton 
proton scattering the two scales are the inverse size of the proton, 
proportional to the QCD A-parameter, and the virtuality, Q, of the 
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photon exchanged between the lepton and the proton. Factorization is 
given by the operator product expansion which separates the hard and 
soft scales while the DGLAP equation governs the Q^— dependence. 

For processes having only one (hard) transverse momentum scale 
the Balitsky, Fadin, Kuraev, Lipatov(BFKL) equation[4,5] determines 
the center of mass energy dependence when the hard transverse mo- 
mentum scale is held fixed. BFKL effects also contribute to two-scale 
processes like deep inelastic scattering, and these effects arc accounted 
for through a resummation of certain higher order terms in the anoma- 
lous dimension and coefficient functions which occur in the operator 
product expansion [6- 12]. 

However, even for processes having only one transverse momen- 
tum scale the BFKL equation ultimately breaks down, at very high 
energies,because gluons diffuse away from the hard scale and reach 
transverse momenta of size A where perturbation theory no longer 
applies [5, 13]. It can be expected that a similar phenomenon will oc- 
cur when BFKL effects are included in the DGLAP equation through 
resummation. For deep inelastic scattering at moderate values of x 
the leading twist terms in the operator product expansion represent 
the structure functions so long as Q/A >> 1. However, if Q/A is fixed 
at some large value the operator expansion breaks down for sufficiently 
small values of x because of diffusion, in both the anomalous dimen- 
sion and coefficient functions, into regions of transverse momentum 
on the order of A. This breakdown of the operator product expansion 
is not due to higher twist terms, but rather to an inability to properly 
separate hard and soft scales at very small values of x. Of course 
within a region of x, say xq < x < 1, one can always choose a Q{xo) 
such that if Q > Q{xo) then the operator product expansion applies. 
As we shall see a little later on (See (40)) this region is determined by 
> C-^^^in x/x]i/3, with b = when a resumma- 

tion of leading logarithmic terms is taken into account, x is a fixed 
parameter, x ~ 1/10 perhaps, not determined in our analysis. 

Our model of diffusion into the infrared is the fixed-coupling BFKL 
equation. Because the fixed-coupling BFKL equation underestimates 
diffusion into the infrared, and overestimates diffusion into the ultra- 
violet, our result on the region of applicability of the operator product 
expansion, (40), should be taken with some care. This uncertainty is 
implicit in (40) in that we are unable to specify the scale of the run- 
ning coupling in D. In the limit for £n ^^'^"^ given just above we have 
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taken a = a{Q{xo)) in (40). 



2 BFKL evolution and its limitations 



The BFKL equation naturally applies to high energy cross sections 
were there is only a single transverse momentum scale, Q, if that scale 
is hard enough so that the running coupling, a{Q), is small. While 
there are a number of processes where BFKL evolution can be directly 
measured, for our purposes the heavy onium-heavy onium total cross 
section is the simplest to consider. The scale Q is determined by the 
inverse heavy onium radius. The cross section is given [14], in leading 
logarithmic approximation, by 



a{Y) = J d^x J\z J (fix j\z'^{x',z')f{x',Y,x)^{x,z) (1) 

where ^{x,z) is the square of the light-cone wave function of the onium 
with X the transverse coordinate separation of the heavy quark and 
heavy antiquark while z is the longitudinal momentum fraction of the 
heavy quark. Y = in s/M"^ with s the center of mass energy squared 
and M the mass of the onium. The BFKL impact parameter ampli- 
tude, f{x', Y, x), is given in terms of the momentum space amplitude 

by 

/(^'' ^' ^) = / 1^ ^ e^-^-^^'-^'f{k', Y, k). (2) 
Asymptotically, 

g(ap-i)y lr?x' Ix 

f{x',Y,x) ~ Anxx'a =exp{ } (3) 

^laN^Cm ^DY 

and 

f(k ,Y,k) ~ — — , = exp{ TTrd— I (4) 

where ap - 1 = ^^In 2 and L> = Using (3) in (1) gives 
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g(ap-l)y 

a - 167ri^^a^— p (5) 
^|aiVeC(3)l^ 

with R the radius of the heavy onium. 

Eq.(5) breaks down at large Y for two separate reasons, (i) Even 
though the heavy oniun state has a smah radius, R, the gluons re- 
sponsible for the growth of (3) and (4) with increasing Y diffuse to 
larger distances, determined by the final exponential factors in (3) 
and (4). When the diffusion reaches momenta as small as the QCD 
A-parameter the whole perturbative approach breaks down, (ii) When 
Y becomes large the cross section clearly grows faster than unitarity 
allows so that new corrections must become important which modify 
simple BFKL behavior. 

Consider first the question of diffusion. Eqs.3 and 4 have been 
derived neglecting the running of the QCD coupling. Let fco = \/R. 
Then the minimum momentum to which gluons diffuse in an onium- 
onium collision in km given by 

ir?ko/km = DY (6) 

as determined from (4) taking [exp{— ^\]^"y)2 ~ ^/^ because dif- 
fusion is maximum at rapidities midway between the two onia. (We 
shall give a more complete derivation of (6) later on.) Then the con- 
dition for the running of the coupling to be negligible over the region 
where gluon momenta diffuse is 

ajkm) - ajko) ^ 
a{km) 



which leads to 



Y << ^ (R\ 
147V,C(3)62 a^{k^) ^ ' 

with h = ^^^l2-K^^ ^'^^ where we have used a{k) = jj^j-py;^- 

If one is willing to use a running coupling in the BFKL equation, 

rapidity values 

Y < - — (9) 
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can be reached. One cannot go beyond the hmit set in (9) without 
having fc^-values less than A. 

The Hmit (8), or (9), is not too bad because unitarity constraints 
are expected to be reached, from (5), when rapidities become of size 

Y < in^. (10) 

For a{ko) very small unitarity constraints become important long be- 
fore diffusion constraints. Thus questions of unitarity can in principle 
be studied completely within the fixed coupling approximation[15] for 
sufficiently heavy onia. However, when a{ko) is of moderate size the 
diffusion limit (9) may be reached before unitarity corrections are im- 
portant. In our leading logarithm discussion the limits (9) and (10) 
are comparable at a(fco) is of size 1/3 to 1/2 although numerical stud- 
ies suggest that one in fact must go significantly beyond (10) to see 
strong unitarity corrections[15]. Thus in most practical circumstances 
diffusion constraints will be reached before unitarity constraints. 



3 Incorporating BFKL evolution in a 
DGLAP formalism 

The asymptotic behavior (4) follows directly from the representa- 
tion[4,5,16] 

T(F, Q/m) = I ^exp{^x(A)y + Xin QVm'} (H) 

where 

r(y,Q//x) = ^/(Q,F,/x), (12) 

with 

X(A) = V'(1)-^V'(A)-^V'(1-A) (13) 

and where the path of integration in (11) goes from A = Aq — i oo 
to Ao + i oo with < Aq < 1. In order to cast this into a form 
corresponding to the operator product expansion it is necessary to 
pick out the leading twist part of (11). If we take Q/ijl> 1 then higher 
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powers of (/x/Q)^ can discarded simply by changing the A-integration 
in (11) to a circle, of radius less than 1, about the origin. Call this 
contour C. Then the leading twist part of T, denoted by L, is 



L{Y,Q/ix) ^ I exp{^x(A)F + A^ngVM'} 

JC ^TT? TT 

= / ^Cnexpi^n^n Q-'lix^ + (n - 1)^} (14) 

and determines the anomalous dimension and coefficient functions [17]. 
The n-integration runs parallel to the imaginary axis and to the right 
of the point n=l. (The anomalous dimension, 7„, is unique in the 
leading logarithmic BFKL approximation. The coefficient function, 
Cn, is not unique but the choice above, where matrix elements are 
normalized to 1 at is natural and convenient for our purposes.) We 
note that when Y is large subleading twist terms are also subleading 
in Y since A = 1/2 is the largest of all the saddle points of x(A). 
Using (14) it is straightforward to find[16] 

_i 7r(n-l) 

and 

where is the function inverse to x- If one writes 

T.-!.™!^)" (17) 



and 

oo 

/^n 7r(n - 1) n-1 



yi V y ^ -^^C \ AT 1 

TT 

7V=0 

then (16) gives 

q(n) =N^iN)_ Using (14) one finds 



N-l . ,fc 

klik + 1)! (iv - 1 - ky. d\k 



N-l-k 



(19) 

A=o 
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with C(0) = 1 and 



p(A) = 2x(A)-i (20) 

(Eq.l9 is most easily found by setting Q = in (14), and then taking 
the term on each side of that equation.) 

The coefficient function and anomalous dimension have a singular- 
ity at n = ap corresponding to the saddle point of xW at A = 1/2. 
One can determine the behavior of Cn and 7^ near n = ap either from 
(15) and (16) using 

x(A) = 2^n2 + 7C(3)(A-^)2 + ... (21) 
or from (17) and (18) using 



(Ain 2)^ / in 2 
nZoo VN Vl47rC(3)- ^ ' 

Near n = ap one finds[4,5,6-8,17] 

'^^ ~ A /n/ ^^^^ 
A^/D{n — ap) 

and 



1 1 (n — ap) 
^"-2"-2V^^- ^''^ 
We are now in a position to see more clearly how diffusion con- 
tributes to L. It is convenient to separate the non-diffusion parts of L 
according to 

L = ^L(g//i) e("--i)^ (25) 
where we have normalized L to be a probability distribution. Then, 



If dn 1 
L{Y,Q/fi)= -J — Cnexp{{jn-^)inQyfi^ +(^n-ap)Y} (26) 

obeys, at large Y, the diffusion equation[4,5,13] 
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as is easily seen using (24) in (26). The solution is 

fa' Q 

^1 ADY 



L(r,Q/M) = '^^V ^ (28) 



which, using (25) and (12), gives (4), 

4 Limitations on the use of DGLAP 
evolution at small x 

Our normal picture of DGLAP evolution is not one of diffusion but 

of a monotonic increase of in Q/^ with Y. In order to see where the 
diffusion is hidden write L, given in (14), as 

MY, Q/^i) = r dyAiy, Q/^)C{Y - y) (29) 
Jo 

where asymptotically. 



VmiDY 



and 



A{Y,Q/^i) = J ^e^"'" QVm'+{"-i)^ 

II Y ^/IGtvDy''^^^ 4.DY ^ ^ ' 

A straightforward calculation shows that the values of y in (29) which 
dominate the integral are 

Y 



^ ^in^ Qlix> 

Thus, for Y » ^" 4^^^ the coefficient function covers most of the 
rapidity region in L, as given in (29), so that for very large values of 
Y most of the diffusion is hidden in C. 
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However, this is not the whole story. There is also significant 
diffusion in A despite the fact that the DGLAP equation 



?^IM.l\,,iy)AiY-y,QM. ,31, 

seems to suggest increasing Q-values and increasing Y-values go to- 
gether. To see where the diffusion is hidden it is convenient to write 

A = 9.e('^p-^')^A{Y, Q/fi) (32) 



So that A obeys 

dA 



din Q/fi 



r dyi^-l/2){y)A(Y-y,Q/ti) (23) 
Jo 



with 



It is now straightforward to see that the y-values dominating (33) are 
given by 

ADY , , 

^ " In^ Q/,. + ADY ^- 

Thus when Y > the action of the anomalous dimension is 

very nonlocal in the DGLAP equation (33). This allows the diffusion, 
which according to (29) is certainly contained in A, to be hidden in 
the anomalous dimension 7. 

To see the general limitations on what values of Y allow a consis- 
tent operator product expansion, and a DGLAP formalism, we note 
that 

L{Y, Q/f,) = -^LiY - y, Q/k)L{y, k/f,). (36) 

Prom (36) one determines , using (25) and (28), that the smallest 
important values of k in (36) are reached when one takes 

Y tnQ/jx 
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with the effective minimum value, km, given by 

en'' ^ = 4DY. (38) 

Taking fe^ = A as the mimimum allowed value for km before the whole 
formalism breaks down we get the restriction 

If in \J? I « in Q/ II one finds Y < j^in'^ Q/n and the breakdown 
occurs because of diffusion in the high order corrections to the anoma- 
lous dimension. U in Q / fi « in ;U^/A^ one finds Y < jjjin'^ ^ I 
and the breakdown occurs in the coefficient function or, in more gen- 
eral schemes, in a combination of the coefficient function and the op- 
erator matrix element. In applying (39) to deep inelastic scattering 
one should identify /v, with the transverse momentum scale where per- 
turbation starts to apply while Y = in x/x where x as the x- value 
where soft gluon emission starts to become important. As a rough 
estimate we take = A and we suppose x 1/10. Then 

in x/x < ■^in'^Q/A (40) 

gives the region of x where the operator product expansion applies. 
(At our level of discussion we are unable to fix the scale at which a 
should be evaluated in D. a = a{Q) would certainly give an upper 
bound in (40).) 

If one wishes to use the DGLAP equation to evolve from Qo to Q 
one must satisfy (39) with /i = Qq. However, the stronger constraint 
is that the operator product expansion be valid at Qq, that is one 
needs in x/x < j^in'^ Qo/A for x- values which are important at Qo 
in determining the parton distribution at the desired small value of x. 
In case one chooses Qo > Q then the strongest constraint is (40), the 
requirement that the operator product expansion apply at the x and 
Q' value one is interested in. 

At the level we have calculated no breakdown of the operator prod- 
uct expansion is yet visible[10,ll]. However, once momenta as low as 
A occur, and that will happen when (39) is not satisfied, higher order 
running coupling corrections will not be small. This signifies the fact 
that the anomalous dimension and/or the coefficient function cannot 
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be calculated perturbatively. Since this breakdown of perturbative oc- 
curs at the leading twist level it corresponds to the breakdown of the 
operator product expansion itself. Also, since our discussion concerns 
leading twist effects this phenomena is different from the appearance 
of renomalons. The breakdown of the perturbation series when higher 
order running coupling effects are included will here not be associ- 
ated with an n! growth of the perturbation series, but simply with 
the fact that higher order running coupling corrections are not small. 
It is likely that diffusion puts significant restrictions on the use of a 
DGLAP analysis of present small-x and small-Q^ deep inelastic scat- 
tering data. At the same time we are presented with the interesting 
challenge of how to analyze moderate and very low x phenomena. 
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